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ABSTRACT 

We have discovered two recycled pulsars in relativistic orbits as part of the first high-frequency survey of inter- 
mediate Galactic latitudes. PSR Jl 157-51 12 is a 44 ms pulsar and the first recycled pulsar with an ultra-massive 
(M > 1.14 M©) white dwarf companion. Millisecond pulsar J1757-5322 is a relativistic circular-orbit system 
which will coalesce due to the emission of gravitational radiation in less than 9.5 Gyr. Of the ~40 known circular 
orbit pulsars, J 1757-5322 and J 1 157-51 12 have the highest projected orbital velocities. There are now three local 
neutron-star/white-dwarf binaries that will coalesce in less than a Hubble time, implying a large coalescence rate 
for these objects in the local Universe. Systems such as Jl 141-6545 (Kaspi et al. 2000) are potential gamma- 
ray burst progenitors and dominate the coalescence rate, whilst lighter systems make excellent progenitors of 
millisecond pulsars with planetary or ultra-low mass companions. 

Subject headings: binaries: close — gammarays: bursts — pulsars: individual (Jl 157-51 12, J1757-5322) — 
stars: planetary systems — stars: white dwarfs — surveys 



1. INTRODUCTION 

The first surveys of the Galactic plane at high frequencies for 
radio pulsars (Clifton & Lyne 1986; Johnston et al. 1992) un- 
covered a much younger and more distant population than the 
traditional 400 MHz surveys. High frequency surveys are much 
less affected by the deleterious effects of scatter-broadening and 
dispersion smearing which make short period pulsars impossi- 
ble to detect at large dispersion measures. Unfortunately, the 
sampling rates employed in these early surveys limited their 
sensitivity to millisecond pulsars (MSPs) (Johnston & Bailes 
1991), and the spatial coverage was limited by the small beam 
of a radio telescope at high frequencies. No millisecond pul- 
sars were discovered in these early high-frequency surveys. It 
was widely believed (Foster, Fairhead, & Backer 1991) that 
MSPs were steep-spectrum objects, and that there was little 
point in searching for them at high frequencies. As the pop- 
ulation of MSPs became larger however, a series of spectral 
studies (Lorimer et al. 1995; Toscano et al. 1998) provided ev- 
idence that spectra of millisecond pulsars were not too dissimi- 
lar from that of normal pulsars. In a simulation of the Galactic 
millisecond pulsar population, Toscano et al. (1998) predicted 
that high frequency surveys undertaken with the Parkes multi- 
beam system (Staveley-Smith et al. 1996) should detect a large 
number of millisecond pulsars. Motivated by this study we have 
conducted the first large-scale survey for millisecond pulsars at 
intermediate Galactic latitudes. 

Our survey covered the region enclosed by 5° < \b\ < 15° 
and -100° < / < 50°. Pointed observations of 265 seconds 
were made with the 64-m Parkes radiotelescope using the 21- 
cm multibeam receiver. The receiver consists of thirteen dual 
linear polarization feeds arranged in a sparse focal-plane pat- 
tern which allows full sky coverage with beams overlapping at 
the half -power points. This increases the rate of sky coverage 
by a factor of thirteen over a single-beam receiver and enabled 
us to conduct a large-scale survey in only 14 days of integration 
time. The system is sensitive to frequencies from 1230 MHz to 
1530 MHz, with an average total system noise temperature of 
24 K. The backend system was built at Jodrell Bank for the 



ongoing deep Galactic plane multibeam pulsar survey (Lyne et 
al. 2000). It includes twenty-six filterbanks, each with ninety- 
six channels and total bandwidth of 288 MHz, centred at a sky 
frequency of 1374 MHz. Filterbank outputs were summed in 
polarisation pairs and one-bit digitised with an integration time 
of 125 (is. The large bandwidth, narrow channels, fast sam- 
pling, low system temperature and multiple beams facilitiated 
the discovery of a large number of pulsars, including millisec- 
ond pulsars, with greater time efficiency than any previous large 
survey. Data were recorded on magnetic tape for offline pro- 
cessing on the Swinburne supercluster, a network of 64 Compaq 
Alpha workstations, using standard techniques (Manchester et 
al. 1996). 

The survey is now complete. Full details of the survey will be 
provided in a subsequent paper (Edwards et al., in prep.). In all, 
we have discovered 58 new pulsars, of which one has a period 
of 44 ms and seven others have periods less than 20 millisec- 
onds. In this Letter we will discuss the two recycled pulsars for 
which we have phase-connected solutions, PSRs Jl 157-51 12 
andJ1757-5322. 

2. PSR Jl 157-51 12 — ULTRAMASSIVE COMPANION 

PSR J 1 157-51 12 is a 44 ms pulsar in a 3.5 d circular orbit, 
discovered in an observation made on 1999 Jan 19th. Its param- 
eters, derived from pulse time-of-arrival analysis, are listed in 
Table 1. The pulsar's relatively short pulse period, small spin- 
down rate and circular orbit indicate recycling in a binary from 
an evolved companion (van den Heuvel 1984). If the compan- 
ion formed a neutron star after recycling the pulsar, we would 
expect an eccentric orbit due to the sudden mass loss associ- 
ated with the neutron star's production. The circularity of the 
orbit implies that the companion is a white dwarf. Since ob- 
served neutron star masses are consistent with a Gaussian dis- 
tribution with a standard deviation of just 0.04 M Q around the 
mean of 1.35 M Q (Thorsett & Chakrabarty 1999), a conser- 
vative lower limit on the mass of the companion white dwarf 
can be obtained by assuming a 1.27 M© pulsar and an edge- 
on orbit. This means that Mwd > 1-14 M©, indicating an 
"ultra-massive" ONeMg degenerate, and could easily be near 
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the Chandrasekhar limit if the orbit is near the median inclina- 
tion angle of binaries. By similar reasoning, if the white dwarf 
mass is less than 1.35 M , the pulsar mass is at most 1.75 M©. 

A number of ultra-massive white dwarf systems are now 
known. The largest sample are those that have been discovered 
in recent X-ray and extreme ultraviolet observations (Bergeron, 
Saffer, & Liebert 1992; Vennes et al. 1997), the masses of 
which are estimated from model-dependent spectroscopic mea- 
surements. The companion to PSR Jl 157-5 1 12 is the first ultra- 
massive white dwarf to have a firm lower mass limit which re- 
lies only on the narrow band of observed neutron star masses 
and Newtonian physics. Steller mergers have been suggested 
for the formation of ultra-massive white dwarfs, however it is 
clear from the maintenance of a binary association that this is 
not the case with the companion of PSR Jl 157-51 12, which 
must have resulted from largely standard stellar evolution. 

PSR J 1 157-5 1 1 2 is the latest example of a pulsar system that 
possesses a very massive white dwarf companion. Recent opti- 
cal observations of the eccentric binary pulsar B2303+46 have 
indicated that a faint object consistent with a white dwarf is co- 
incident with the position of the pulsar (van Kerkwijk & Kulka- 
rni 1999). If this proves to be a genuine association, the peri- 
astron advance and mass function limit the companion mass 
to greater than 1.2 M Q . The newly discovered pulsar binary 
Jl 141-6545 (Kaspi et al. 2000), with a minimum companion 
mass of 0.95 M Q is probably also a white dwarf (Kaspi et al. 
2000; Tauris & Sennels 2000). The orbits of these systems are 
highly eccentric and thus it is believed that the neutron stars 
were born after their white dwarf companions. Since the pulsar 
in the PSR J 1 157-51 12 binary is recycled, it is the first exam- 
ple of an ultra-massive white dwarf that was formed after the 
neutron star. 

The high mass ratio between any possible progenitor to the 
white dwarf and the pulsar must have resulted in a phase of 
unstable mass transfer. The system probably passed through a 
common envelope phase when the white dwarf progenitor was 
on the asymptotic giant branch (AGB), as suggested by van den 
Heuvel (1994) for the PSR J2145-0750 system. Recent mod- 
els (Ritossa, Garcfa-Berro, & Iben 1999) indicate that main se- 
quence stars as massive as 1 1 M Q can avoid core collapse on the 
AGB to produce ultra-massive white dwarfs (Mwd = 1-37 M Q 
in the above case). Given that envelope loss in the common en- 
velope phase reduces the resultant remnant mass from a donor 
star of given mass (Iben & Tutukov 1993), we expect that the 
progenitor to the white dwarf companion of PSR Jl 157-51 12 
weighed in excess of 10 M Q . This would mean that the com- 
panion's radius (Ritossa, Garcfa-Berro, & Iben 1996; Garcfa- 
Berro, Ritossa, & Iben 1997; Iben, Ritossa, & Garcfa-Berro 
1997; Ritossa, Garcfa-Berro, & Iben 1999) was < 800 R Q at 
the time of Roche lobe overflow. The orbital energy lost during 
spiral-in to the current orbit was insufficient to expel the en- 
velope of the companion, without an additional energy source 
such as an accreting neutron star. Future proper motion mea- 
surements will help constrain models for its evolution. For an 
11 M Q companion with a pre-CE Roche lobe radius > 500 
R Q (requiring that at most 7/8 of the envelope ejection energy 
comes from non-orbital sources), the kick velocity of the pulsar 
birth event must have been < 120 km s" 1 for the system to have 
remained bound. On the other hand, if the system was much 
closer at the time of explosion much greater kicks could be ac- 
commodated and we might expect the system to have a very 
large runaway velocity (Bailes 1989). 



Along with J1757-5322 (below), Jl 157-51 12 brings the to- 
tal number of intermediate mass binary pulsars (those with 
white dwarf companions heavier than 0.6 M Q ) (Arzoumanian, 
Cordes, & Wasserman 1999; Tauris & Savonije 1999; Camilo 
et al. 1996) to seven. The companion of PSR J 1 157-51 12 is 
by far the most massive of the class. That the neutron star sur- 
vived what seems to have been a common envelope phase with 
a massive AGB star strongly supports the case against black 
hole formation via hypercritical accretion (Armitage & Livio 
2000). 

3. PSR J1757-5322 - COALESCING BINARY 

PSR J1757-5322 is an 8.9 ms pulsar which was discovered 
in an observation made on 1999 May 8th, at an epoch where 
the orbital acceleration experienced by the pulsar (~7.5 km s" 2 
in the line-of-sight) resulted in noticable period evolution in the 
265-s integration. Subsequent observations revealed the pulsar 
to be in a nearly circular 1 1-h orbit with a companion of at least 
0.55 M Q . A phase coherent timing solution has been obtained, 
the parameters of which are listed in Table 1 . From similar rea- 
soning to that invoked above for Jl 157-51 12, we expect that 
the companion is a white dwarf. Its parameters are compati- 
ble with the common envelope model discussed above for the 
intermediate mass binary pulsars. 

The J1757-5322 system is highly relativistic, having the 
highest orbital velocity of all circular orbit binaries (followed 
by Jl 157-51 12). The near-circularity of the orbit reduces the 
measurability of relativistic effects, however some effects could 
be measurable in the coming decades with long timing base- 
lines. Specifically, assuming an inclination angle ;' = 60° (the 
median value for randomly oriented orbits), General Relativ- 
ity predicts that the periastron of J1757-5322 is advancing at 
a rate of ~ l.°2 yr -1 , that its orbital period is decreasing by 
1.8 /xs yr" 1 , that the pulsar's spin axis is precessing an entire 
turn every 1700 yr and that the system will coalesce in at most 
9.5 Gyr due to the emission of gravitational waves. The Shapiro 
delays due to space-time curvature around the companion (for 
; = 60°) are 29 [is and 13 [is peak-to-peak respectively for PSR 
Jl 157-51 12 and PSR J1757-5322, and may be considerably 
longer if the systems are nearly edge-on to the line of sight. 

After about 7.6 Gyr the orbital period of J1757-5322 will 
have decreased to one-third its current value, resulting in four 
times the orbital acceleration. Such a system, if observed in our 
survey would have experienced significant signal-to-noise ratio 
loss, limiting the chance of detection. However, such a system 
will only live for ^400 Myr, as compared to ^8 Gyr for the 
preceding less accelerated post-CE evolution phase in which 
J1757-5322 was detected. Such highly accelerated systems are 
thus likely to be an order of magnitude lower in detectable pop- 
ulation. 

The companion of J1757-5322 will fill its shrinking Roche 
lobe when it has an orbital period of ^50 s. The range of poten- 
tial outcomes of this is diverse and interesting. Models suggest 
that after a period of mass transfer, the white dwarf will tidally 
disrupt and form an accretion disk around the pulsar in just a 
few seconds if the companion mass is greater than about 0.7 
M Q (van den Heuvel & Bonsema 1984). The rapid mass trans- 
fer will result in an enormous release of energy, potentially in 
the form of a gamma-ray burst. However, it is not clear whether 
or not the accretion luminosity would be self-limiting due to 
the outward radiation pressure inhibiting mass transfer (Fryer 
et al. 1999). If there is no remnant of the white dwarf we might 
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expect a solitary millisecond pulsar with a very rapid rotation 
rate (van den Heuvel & Bonsema 1984) (possibly like PSR 
B 1937+21) and bigger than average mass. Alternatively it may 
be that the mass of the pulsar is driven over the critical limit for 
neutron stars and becomes a black hole. If the disk surrounding 
the neutron star is not completely destroyed the formation of a 
planetary system around a millisecond pulsar is not unreason- 
able, similar to PSR B 1257+12. On the other hand if just 5% 
of the mass of the white dwarf is retained intact, a short-period 
eclipsing pulsar system such as PSR J205 1-0827 may result. 
The measurement of the proper motion of the system is there- 
fore vital in addressing these issues. A large velocity would 
provide a good match to the planet pulsar PSR B 1257+12. A 
very small velocity is probably less conclusive. 

The observed rate of gamma-ray bursts in the local Universe 
(Phinney 1991) is now broadly compatible with the expected 
coalescence rate of pulsar-white dwarf binary systems. There 
are now three such systems known that will coalesce in the 
lifetime of the universe. PSR J0751+1807 will coalesce in 
~7 Gyr. It is only 2 kpc away but has a light 0.15 M Q com- 
panion that upon reaching its critical Roche lobe will undergo 
stable mass transfer to its much heavier pulsar companion, re- 
sulting in an ultra-compact X-ray binary (Ergma, Lundgren, & 
Cordes 1977). PSR Jl 141-6545, as mentioned above, has a 
white dwarf companion of at least 0.95 M Q in an eccentric or- 
bit that will coalesce in only ~1.3 Gyr, and is ^3.2 kpc distant. 
It is a prime candidate for rapid mass transfer after the onset of 
Roche lobe overflow and a potential gamma-ray burst progen- 
itor due to its high mass. The alarming thing about this pulsar 
is its youth. Its characteristic age is just 1 Myr and it may not 
pulse for much more than 10 Myr due to magnetic dipole break- 
ing pushing it past the rotation limit for the emission of radio 
pulses. If we take the standard pulsar beaming fraction of ~5, 
and the short lifetime of the binary in an observable state (~10 
Myr) compared to the coalescence time (1.3 Gyr), the extrap- 
olated Galactic coalescence rate for such objects is enormous! 
We could reasonably expect there to be 5 x 1300/10= 850 dead 
white dwarf neutron star binaries within a few kpc of the Sun 



that will coalesce in less than a Hubble time. The total Galactic 
population would then be 1 -2 orders of magnitude greater than 
this, and the coalescence rate near 10" 5 yr" 1 per Galaxy, similar 
to that expected of double neutron star binaries, another candi- 
date for gamma-ray bursts (Phinney 1 99 1 ; Narayan, Paczyhski, 
&Piran 1992). 

The contribution of PSR J1757-5322 to the inferred 
deathrate of white dwarf-neutron star binaries is much smaller. 
With the relative proximity (~2 kpc) but enormous character- 
istic age (~5 Gyr) of PSR J1757-5322, we might expect ~250 
similar coalescing systems in the Galaxy, with a coalescence 
rate of 5 x 10" 8 yr _1 if we assume an evenly distributed disk pop- 
ulation 15 kpc in radius. Hence there should be ^500 systems 
in the Galaxy which have already coalesced, compatible with 
the observed population of Galactic disk pulsars with planetary 
or evaporating white dwarf remnant companions (three within 

1.5 kpc). 

The white dwarf companions of PSR J 1 157-51 12 and PSR 
J1757-5322 are possibly observable with 10m class telescopes 
or HST PSR B0655+64 (M WD > 0.64 M Q ) has a timing age of 

3.6 Gyr and a 22nd R-magnitude companion (Kulkarni 1986) 
at a distance of about 500 pc. The pulsars described here are 
probably 2-4 times as distant and of comparable timing ages 
so we might expect their companions to have R-magnitudes of 
24-26. Detections of the companions to PSR J1757-5322 and 
PSR J 1 157-51 12 will provide invaluable information about the 
cooling times of massive white dwarfs as the pulsars give an 
independent estimate of their age. 
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an APA scholarship (R.E.) and ARC fellowships (M.B.) and 
grants. 



REFERENCES 



Armitage, P. & Livio, M. 2000, ApJ, 532, 540. 

Arzoumanian, Z., Cordes, J. M., & Wasserman, I. 1999, ApJ, 520, 696. 
Bailes, M. 1989, ApJ, 342, 917. 

Bergeron, P., Saffer, R. A., & Liebert, J. 1992, ApJ, 394, 228. 

Camilo, E, Nice, D. J., Shrauner, J. A., & Taylor, J. H. 1996, ApJ, 469, 819. 

Clifton, T. R. & Lyne, A. G. 1986, Nature, 320, 43. 

Ergma, E., Lundgren, S. C, & Cordes, J. M. 1977, ApJ, 475, L29. 

Foster, R. S., Fairhead, L., & Backer, D. C. 1991, ApJ, 378, 687. 

Fryer, C, Woosley, S., Herant, M., & Davies, M. 1999, ApJ, 520, 650. 

Garci'a-Berro, E., Ritossa, C, & Iben, I., J. 1997, ApJ, 485, 765. 

Iben, I. & Tutukov, A. 1993, ApJ, 418, 343. 

Iben, I., J., Ritossa, C, & Garci'a-Berro, E. 1997, ApJ, 489, 772. 

Johnston, S. & Bailes, M. 1991, MNRAS, 252, 277. 

Johnston, S., Lyne, A. G., Manchester, R. N., Kniffen, D. A., D'Amico, N., 

Lim, J., & Ashworth, M. 1992, MNRAS, 255, 401. 
Kaspi, V. M. et al. 2000, ApJ, . in press, astro-ph/0005214. 
Kulkarni, S. R. 1986, ApJ, 306, L85. 



Lorimer, D. R., Yates, J. A., Lyne, A. G, & Gould, D. M. 1995, MNRAS, 
273,411. 

Lyne, A. G. et al. 2000, MNRAS, 312, 698. 
Manchester, R. N. et al. 1996, MNRAS, 279, 1235. 
Narayan, R., Paczyhski, B., & Piran, T. 1992, ApJ, 395, L83. 
Phinney, E. S. 1991, ApJ, 380, L17. 

Ritossa, C, Garci'a-Berro, E., & Iben, I., J. 1996, ApJ, 460, 489. 

Ritossa, C, Garci'a-Berro, E., & Iben, I., J. 1999, ApJ, 515, 381. 

Staveley-Smith, L. et al. 1996, Proc. Astron. Soc. Austral., 13, 243. 

Tauris, T. M. & Savonije, G. J. 1999, A&A, 350, 928. 

Tauris, T. M. & Sennels, T. 2000, A&A, 355, 236. 

Thorsett, S. E. & Chakrabarty, D. 1999, ApJ, 512, 288. 

Toscano, M., Bailes, M., Manchester, R., & Sandhu, J. 1998, ApJ, 506, 863. 

van den Heuvel, E. P. J. 1984, J. Astrophys. Astron., 5, 209. 

van den Heuvel, E. P. J. 1994, A&A, 291, L39. 

van den Heuvel, E. P. S. & Bonsema, P. T. J. 1984, A&A, 139, L16. 

van Kerkwijk, M. & Kulkarni, S. R. 1999, ApJ, 516, L25. 

Vennes, S., Thejll, P. A., Galvan, R. G, & Dupuis, J. 1997, ApJ, 480, 714. 



Table 1 

astrometric, spin, binary and derived parameters 





Jl 157-51 14 


J1756-5322 


Right ascenscion, a (J2000.0) 


ll h 57 m 08. s 166(1) 


17 h 57 m 15?1618(4) 


Declination, S (J2000.0) 


-51°12'56."14(3) 


-53°22'26"38(1) 


Pulse period, P (ms) 


43.58922706284(12) 


8.869961227275(4) 


P epoch (MJD) 


51400.000000 


51570.000000 


Period derivative, P (10~ 20 ) 


14.6(16) 


2.78(15) 


Dispersion Measure (pc cm" 3 ) 


39.67(3) 


30.793(4) 


Orbital Period, P orb (d) 


3.50738639(3) 


0.4533112382(7) 


Projected semi-major axis, a sin; (lt-s) . . 


14.28634(3) 


2.086526(5) 


Epoch of Ascending Node, r asc (MJD) . 


51216.4442642(14) 


51394.1080692(3) 


ecosw a 


-0.000322(4) 


-1.0(40)xl0- 6 


esinw 


0.000240(4) 


4.3(45) xlO" 6 


Derived Parameters 


Longitude of periastron, uj (°) 


306.7(6) 


347(58) 


Orbital eccentricity, e 


4.02(4) xlO" 4 


4.4(45) xlO" 6 


Companion mass, Mwd (M q ) 


>1.14 


>0.55 



Characteristic age, t c (Gyr) 4.70 5.04 

Surface magnetic field, B surf (10 8 Gauss) 25.55 5.02 



a We used the ELL1 binary timing model of TEMPO (Wex, unpublished). To avoid the 
large covariance between the time of periastron (7b) and w for e <C 1, the time of ascending 
node (r asc ) and the Laplace parameters ecosw and esinw are used instead. 

Note. — Values in parentheses apply to the final digit of the value and represent twice 
the formal uncertainties produced by TEMPO after scaling TOA uncertainties to achieve a 
reduced x 2 of unity. 



